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Ziegler et al. reply: In their Comment jjj on our re- 
cent Letter, Nersesyan and Tsvelik (NT) have raised 
questions about the applicability of our calculation of the 
density of states (DOS) of a 2D d-wave superconductor 
(SC), though they do not actually dispute the result of a 
nonzero Fermi Level DOS for a d-wave SC. They question 
the relevance of our exact calculation of the DOS using a 
Lorentzian disorder distribution, claiming "that models 
with Lorentzian and Gaussian disorder belong to differ- 
ent universality classes" , on grounds of a comparison with 
a straightforward perturbation expansion in the disorder 
potential. We use a different method ^,|| to obtain lower 
bounds for the DOS for other (e.g. Gaussian) distribu- 
tions supporting our contention of a nonzero DOS at the 
Fermi level. NT challenge this method by alleging that 
"similar discrepancies exist between the results obtained 
by this method for another model of disorder - (2+l)D 
fermions with random mass" , i.e. they contend that the 
method used in (J] gives incorrect results. 

We have prepared a paper pj which explicitly proves 
the existence of a nonzero lower bound for the Fermi 
level DOS for 2D SC's with nodes (like p-wave, d-wave 
or extended s-wave SC's) for a large class of disorder 
distributions, including the Gaussian. The proof applies 
the work of Ref. || to the case of interest. It also shows 
that the contribution of the "tails" of the distribution 
are unimportant for the Fermi level DOS. Therefore, our 
method of dealing with disorder in our Letter (using a 
Lorentzian distribution) produces generic results for the 
physical system of a 2D SC with nodes, with disorder 
in the chemical potential. NT have not identified an er- 
ror in our proof. Their argument against it is based on 
inconclusive numerical evidence, as discussed below. 

The approach of NT |^| approximates a lattice sys- 
tem with disorder in the chemical potential by a contin- 
uum model of Dirac fermions with random gauge fields. 
The more general model of Dirac fermions with all kinds 
of disorder allowed (e.g. random mass) has been ana- 
lyzed by Mudry et al. Rj who found that the critical 
line corresponding to pure gauge field disorder is unsta- 
ble with respect to the other kinds of disorder. They 
conclude that "unless there exists a symmetry of the un- 
derlying lattice regularization" forbidding other disorder 
than random gauge fields the non-Gaussian cumulants 
are relevant perturbations. In short, the model consid- 
ered by NT in Ref. || is most likely not generic, since 
the physical system does not provide a symmetry that al- 
lows the disorder to be described solely by random gauge 
fields. In this regard we would like to emphasize that our 
work H is not "aimed to debunk the entire field theoret- 
ical approach to disordered systems", but is actually in 
agreement with the conclusions of the no less field theo- 
retical work by Mudry et al. M . 

1) Comparison with perturbation theory. It is 

clear that the Lorentzian disorder distribution for which 
we obtained our exact results is quite special in terms of 


perturbation theory, but we do not believe that for the 
Fermi level DOS it will give qualitatively different results 
than the impurity model with Gaussian disorder consid- 
ered by NT. The fact that the Lorentzian disorder leaves 
the DOS of an s-wave SC unaffected shows that its tails 
do not necessarily lead to unphysical results. Explicit 
calculations show that the tails of the distribution 
are in fact unimportant if the DOS at the Fermi level 
is considered. In any case, it is disingenuous to argue 
that because we can not reproduce a perturbative cal- 
culation with a Lorentzian distribution the result of our 
Letter must be "in a different universality class". Since 
every moment of the distribution is infinite, perturbative 
arguments simply do not apply to this case. 

2) Dirac fermions with random mass. In contrast 
to what NT imply, there is analytical ||,|| and numeri- 
cal J8|,|9| evidence that Dirac fermions with random mass 
have a nonzero DOS at the Fermi level. For a Gaussian 
disorder distribution with width 7 the DOS at the Fermi 
level was estimated in to be oc exp(— tt/j), a nonana- 
lytic behavior with respect to the disorder that can not 
be obtained within perturbation theory. The related 2D 
random bond Ising model is a highly controversial field 
of research JTQ3 . Different groups have obtained different 
results, e.g. Refs. 11 and 12 of the Comment by NT 
disagree with their Ref. 8 on almost everything but the 
specific heat. The numerical work quoted by NT as Ref. 
10 agrees partially with their Refs. 8, 11 and 12 as well 
as with the result of Ziegler j|. Because of the extreme 
weakness of the log log(T — T c ) divergence the numerical 
data are not able to resolve the question of divergence or 
finiteness of the specific heat; the data can be fitted to 
both theories [jfl]]. Therefore, there is no evidence that 
the methods of Refs. produce incorrect results. 
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